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The composite torsional ultrasonic oscillator, a versatile experimental system, can be used to inves- 
tigate slip of Newtonian fluid at a smooth surface. A rigorous analysis of slip-dependent damping 
for the oscillator is presented. Initially, the phenomenon of finite surface slip and the slip length 
are considered for a half-space of Newtonian fluid in contact with a smooth, oscillating solid surface. 
Definitions are revisited and clarified in light of inconsistencies in the literature. We point out that, 
in general oscillating flows, Navier's slip length b is a complex number. An intuitive velocity disconti- 
nuity parameter of unrestricted phase is used to describe the effect of slip on measurement of viscous 
shear damping. The analysis is applied to the composite oscillator and preliminary experimental work 
for a 40 kHz oscillator is presented. The Non-Slip Boundary Condition (NSBC) has been verified for 
a hydrophobic surface in water to within ~ 60 nm of \b\ = nm. Experiments were carried out at 
shear rate amplitudes between 230 and 6800 s _1 , corresponding to linear displacement amplitudes 
between 3.2 and 96 nm. 
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1. INTRODUCTION 



In fluid mechanics, the Non-Slip Boundary Condition (NSBC) states that fluid at the interface 
between a fluid and a solid surface is stationary with respect to that surface. Surface slip arises when 
there is non-zero relative motion between the interfacial fluid and the solid surface, in which case the 
NSBC does not apply. In the vast majority of experiments described historically, the NSBC has been 
an adequate boundary condition Q, Q{. However, finite surface slip has been the subject of much recent 
interest and promises to play^ ^significant role in the current and future development of micro- and 

Navier 3] first described surface slip as early as the 1820s and 



nanofluidic devices 



;es to pL 



lis mathematical definition of the 'slip length' has been widely used, for example in recent reviews 



flay,y,a 



and analyses of slip-dependent flow 
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Slip is a well-studied phenomenon in non-Newtonian fluids such as polymer solutions 
[r] ]. For slip of Newtonian fluids (Newtonian slip), theory is well developed Q|, but interpretation of 
experiments has been much more controversial. Newtonian slip is important because technologies can 
usefully be applied to Newtonian fluids, but also because understanding of Newtonian slip will provide 
a more profound understanding of fundamental slip mechanics. The physical mechanisms governing 
slip are affected by variables that include surface chemistry, shear rate, surface structure, Newtonian 
viscosity, non-Newtonian behaviour at the molecular scale and combinations of each of these variables 



Nanoengineered 'superhydrophobic' surfaces trap gas between the fluid and the surface, producing 
large slip lengths of the order ~10 /im 

Experimentally, Newtonian slip has been investigated using a variety of methods [lj]. The most 
widely used and precise technique involves measurement of the drainage forces between closely-spaced 
surfaces in an Atomic Force Microscope (AFM) or Surface Force Apparatus (SFA) and applying an 



analysis based on Reynolds' theory of lubrication 



20l | . Robust, non-zero values of b for Newtonian 



fluids at relatively smooth surfaces are typically of the order 10-100 nm 
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Experiments have not established firm quantitative correlations between slip and major experimental 



variables, such as contact angle 



HQ 



Recently, it has emerged with some clarity that in AFM or 



SFA experiments, the NSBC can be expected to hold on smooth hydrophilic surfaces 
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while finite slip lengths of order 10-20 nm have been measured on some hydrophobic surfaces 
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Additionally, the long-held idea that surface roughness should inhibit slip [2J, |3l[ is now supported by 



reasonable experimental evidence 



ay, 



Shear-mode oscillations can be used to drive fluid flow, and such systems pose scientific questions that 
are distinct from those relating to the study of steady fluid flows. Oscillating piezoelectric components 
are likely to prove very useful in future applications. Ultrasonic oscillators, well established as important 
tools for investigating thin films and bulk fluids, can also be used to probe the interfacial forces that are 



relevant in the study of slip. Kanazawa and Gordon 



32i | analytically coupled a half-space of Newtonian 



fluid to a piezoelectric Quartz Crystal Microbalance (QCM), in which a so 



shear mode, u sing the NSBC. Further reports have presented experiments 

10, fit 
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id interface oscillates in 
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361 ] and analysis 
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[39[ describing slip of a Newtonian fluid adjacent to an oscillating surface. A relatively 



large set of slip lengths measured using a QCM was obtained by Ellis and Hayward 



3J. 



Torsional quartz crystal oscillators operating at kHz frequencies were initially developed and used for 



viscoelasticity measurements 



401 ] and have been used to measure the complex and viscoelastic behaviour 



of polymer solutions. More recent work has investigated dispersions of colloidal 



41 



42 



431 ] and micellar 



44| spheres. The shear rates accessible using torsional oscillators are considered advantageously high for 



such rheological measurements 



43|, 



44|. A particular configuration, the composite torsional ultrasonic 



oscillator, was developed and analysed by Robinson et al. 



451 ] and Robinson and Smedley [4g], who also 



derived an equivalent circuit. Recently, the first experimental results investigating Newtonian slip at a 
smooth surface using such a device have been reported j^], for oscillations at 40 kHz. 

There has been some confusion relating to fundamental aspects of slip, perhaps best characterised by 
Neto et al. [lj], who noted that "dialogue between the QCM community and other researchers interested 
in slip has been scarce". In order to address this issue, Navier's definition of the slip length is used 
consistently in the current work. An intuitive, consistent and general methodology for slip measurement 
at an oscillating surface is presented. Analytic inconsistencies in the literature are identified and clarified, 
especially those concerning an oscillating interface. The composite torsional oscillator is introduced along 
with a full analysis for measurement of Newtonian slip using this device. Preliminary experiments are 
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presented to demonstrate how this device is used to probe interfacial forces in oscillatory flow. We 
discuss the interpretation of measurements and the analytical assumptions on which they are based. 
Finally, methods of slip measurement are compared and the advantages of using the composite torsional 
oscillator are explained. 

2. THEORY 

A. Newtonian slip at a smooth, oscillating surface 

1. Flow profile for shear mode oscillations 

By definition, viscosity r] is constant and independent of shear rate throughout a Newtonian fluid. The 
velocity field u of an incompressible Newtonian fluid of density p in the presence of pressure gradient 
Vp and body force F can be determined by solving the Navier-Stokes equation, 



p { + u.Vu ) = F - Vp + r?V 2 u, 

V.u = 0. 

Microfluidic flows have characteristically low Reynolds numbers, so the convective derivative pu.Vu is 
insignificant when compared with the viscous term t/V 2 u. Solutions to the Navier-Stokes equation are 
dependent on flow geometry and boundary conditions. 

The relevant solution for the current paper concerns a smooth surface oscillating in shear mode in 
contact with a fluid half-space. In the absence of significant body forces or pressure differentials, the 
equation of the shear wave induced in the fluid is 



u f (z,t) = u f0 (i)exp f -^y z \ > ( 2 ) 

where Uf is the fluid velocity in the direction of the surface oscillations, z is the Cartesian direction 

i 

perpendicular to the solid surface and 5 = (^j^j 2 ■ The value of uj at the interface is Ufo, which 
depends on the fluid-surface boundary condition. The time dependence of uj and Ufo is e tujt , where the 
frequency of oscillation is u> and t is time. This solution is relevant for ultrasonic oscillators. 



4 



2. Navier's slip length 

Boundary slip occurs when the NSBC is contravened, so that there is finite relative motion between a 
solid surface and an adjacent fluid. Navier [3] characterised the slip boundary condition by defining the 
slip length b. For flow adjacent to a planar, impermeable solid surface, 



u\\q (t) = b- 



z=0- 



(3) 



dz 

Here uu is the shear velocity relative to and parallel with the surface and um is the value of un at 
z = 0. Navier's definition of slip has been widely used over a long period and has a strong physical basis 
relating to viscous force at the interface. The fluid shear rate at the interface that appears in Eq. [3] also 
determines the viscous force, 

du\\(z,t) 

Vxz {t) = V 



Slip therefore determines the shear stress at the interface, which is analogous to friction in character. 
The coefficient k, defined as 

&xz (t) = hup (t) , (5) 
is finite only when there is slip. Using Eqs. [3l U] and El we find 

6 . I (6) 

It is noted that the NSBC and boundary slip are, respectively, analogous to the regimes of 'static' 
and 'sliding' solid-solid friction. Extending the analogy, we might expect a wide variety of material 
and interfacial properties to determine the onset and magnitude of solid-liquid friction. If slip length 
depends on shear rate, Eq. [3] becomes non-linear and rate-dependent formulations of slip (e.g. [481 ] ) 
become relevant. Currently, there is conflicting experimental evidence relating to possible shear rate 
dependence of 6 Q, 0] . 
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A different physical interpretation of slip length relates b to the distance into the solid over which 
the tangent to the flow profile at the surface must be extrapolated for the velocity to reach zero. In this 
context, Navier's slip length has correctly been described as a fictional distance 3, 



in Fig. [TJ The definition of slip length used in some recent slip studies 
by Ellis and Hay ward [33J as (their notation) vp(dp) = vsidp — b), or 



33j, 



37|, 



38 



491 and is labelled 
□ 

39J, and formalised 



«||(-6 1} t) = 0, (7) 

is distinct from Navier's definition in Eq. [3l Ellis and Hayward's definition corresponds to extrapolation 
of the flow profile a distance b\ (Fig. [T]) relative to the surface when slip occurs. The value of fr i is equ al to 



flfl 



it 



33, 



37|, 



5QC, 



b only in the case of a linear velocity flow profile. An often- used schematic diagram 
which uses a linear flow profile to depict slip length, appears to have been an important factor in the 
emergence of Eq. [7J 



3. An intuitive slip parameter 

We wish to point out that when the flow of a Newtonian fluid is driven by an oscillating surface (as 
introduced above), the shear rate introduces non-zero phase into the value of b, which is defined using 
■U|| (Eq. [3]). The existence of this phase highlights the limitations of direct quantitative comparison of 
measurements of b using different experimental techniques. 

When measuring slip at an oscillating surface, it is both intuitive and useful to use the dimensionless 
complex parameter a = \a\e l< ^ a . This parameter, which simply quantifies the velocity discontinuity 
between the solid surface and the first liquid layer, simplifies the algebra describing the effect of damping 
on the equivalent circuit of a piezoelectric oscillator. It is defined by interfacial properties and does not 
depend on the flow profile, 



au s (t) = u f0 (t) . 



(8) 



The velocity of the solid surface is u s and the fluid oscillation is assumed to be dominated by a sinu- 
soidal oscillation of angular frequency equal to that of the solid surface. Under this assumption, a is 



6 



351 ] and used elsewhere 34j, |5l( for modelling 



equivalent to parameters introduced by Ferrante et al 
of slip at an oscillating surface; to our knowledge, it has not previously been applied to experimental 
measurements with full consideration of phase dependence. More generally, a can be decomposed int o 
Fourier components, an approach that has been used when modelling nanobubble-induced slip 
The NSBC is consistent with the values \a\ = 1 and <b n = 0. 



52, 



5m- 



For an oscillating surface, flow is relative to surface motion, so 



U]| (t) = (u f (t)-u s (t)) 

/2i X X W 



= u s (t) ^aexp I —zj - 1 

The magnitude and phase of u\\ relative to u s are plotted in Fig. for both slip and non-slip cases. 



Figure 2(a) shows that the relative amplitude of oscillation is uniformly non-zero at the surface when 
slip is present, but varies significantly near the interface (| < 0.25) depending on the values of |a| and 
4> a . Further from the surface (| > 1), the relative amplitude tends towards a stable value close to 1 but 



dependent on \a\. The relative phase follows a similar trend in Fig. 2(b) , showing strong dependence on 
the combination of \a\ and (fi a near the interface. There is a phase discontinuity at | as 0.1 for the cases 
in which \a\ > 1. The plane of discontinuous phase is consistent with a negative slip length (Fig. [3|) and 
could, for example, be interpreted as the edge of a stagnant boundary layer. Further from the surface 
(§ 1)> the relative phase converges towards the NSBC value, which tends towards zero a long way 
from the interface. 

Using this flow profile, slip length can be directly related to a (Eqs. [21 EJ [8] and [9]): 

a = —=^. (10) 



Equation [TU] is plotted in Fig. fusing the notation b = |&|e J< V Figure 3(a) shows that \a\ decreases when 
4>b is within 90° of —45°, so that the slip length b is positive. The phase plot (Fig. |3(b)[ ) demonstrates 
that 4>b = —45° when the fluid oscillates in phase with the surface (4> a = 0°). For the NSBC, b = 0, 
so (fib can take any value. The maximum size of the phase shift (fi a is limited by slip length magnitude. 
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The complicated nature of the flow profile near the interface supports the use of the parameter a, which 
does not depend on the flow other than at the interface itself, in contrast to b. However, calculation of 
the shear rate using the flow profile is still necessary for derivation of a from surface force experiments. 

4- Effect of slip on damping force and mechanical impedance 

The damping force on an oscillating surface is equal and opposite to the viscous force in the adjacent 
Newtonian fluid. When slip occurs, the damping stress is (Eqs. [5J H] and EI]), 

■\/2i 

o xz (t) = -^--qaus (t). (H) 

Analysis of ultrasonic oscillators is greatly simplified by using the equivalent circuit analogy, in 
which mechanical components are replaced by their electronic equivalents. We now incorporate the 
damping stress into the analysis of a mechanical circuit. The stress at the interface (Eq. fTTj) . which has 
components both in phase and out of phase with surface velocity, can be written as 



o xz (t) = -(R + iX) u s (t) . (12) 
The impedance of the system is Z = R+iX, as previously applied to a torsional oscillator by Bergenholtz 



et al. [4l|], who used the NSBC. This impedance is dependent on the liquid, the surface, and the 
interaction between them. The real part of the impedance (R) determines the stress opposing, and in 
phase with, the surface velocity, 

R= V^/TTX (i3) 



5 V4 

This impedance is viscous in nature and dissipative of energy. In an equivalent circuit, velocity is 
analogous to current and R is like a resistance. Energy dissipation is determined by the quality factor 
of the oscillation. 

The inertial term (X) gives the stress that is in phase with the acceleration of the surface, 
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x = v^ (i4) 



5 U 

An increase in the inertial term is equivalent to the addition of mass to the surface of the oscillator. In 
the equivalent circuit, the time derivative of current is analogous to acceleration (modified by a factor 
of lo), so X is an inductance. A change in X is measured by considering the change in the period of 
oscillation. Under the NSBC, R and X are equal for Newtonian fluid damping. Slip-induced change in 
the relative sizes of these terms indicates a finite value of <j> a . 

5. Conceptual issues 

The theory that has been presented implicitly addresses two conceptual pitfalls that have caused con- 
fusion and misinterpretation in a range of previous work. The first pitfall concerns the variant 'slip 
length' b±, defined in Eq. O which is only equal to b in the limit of a linear flow profile. This limit is 
generally not applicable for slip measurements. The difference between b and b\ is greater than 10% 
for typical slip measurements using AFM, SFA or ultrasonic oscillation methods. We further note that 
Ellis and Hay ward's definition has been used inconsistently; both with and without considering a Taylor 



expansion 



38, 



39J. Navier's slip length b should be used at the present time due to its strong physical 
basis and to promote consistency across experimental studies. 

The second conceptual pitfall relates to those experimental studies of slip at oscillating surfaces 

3^], especially the assumption that 



33 



33, 



38 



which have imposed a restriction on the phase of slip 
the slip length b is real-valued. It is equivalent to assume that k or its reciprocal s is real-valued, or 
more generally, that the damping force is in phase with the surface velocity (Eq. [6]). Several lines of 
reasoning lead us to emphasize that b must be allowed to have non-zero phase. We have shown that, 
because the slip length is defined using the shear rate (Eq. [3]), the z-dependence of the fluid velocity (a 
complex exponential, Eq. [2]) introduces phase to the value of b. The usual (NSBC) form of solid-fluid 
friction force is viscous (Eq. [TT]) and therefore 45° out of phase with the velocity above an oscillating 
surface. Figure [3] shows that the phase of b in the limit approaching the NSBC is non-zero and that 
both magnitude and phase must be determined for either b or a to be fully described. Restricting b to 
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real values is equivalent to restricting a to one degree of freedom in Eq. ITOl Furthermore, experiments 
[47I have indicated that a fluid can be expected to oscillate in phase with, or slightly lag, the surface 
oscillation. Neither interpretation is consistent with a real value of b. 

Having established that b must be considered complex, we address specific arguments made in pre- 
vious studies in support of a real-valued b. Firstly, there has been a direct assumption 



38 



39] 



that the value of k (or s) should not be phase dependent. This assumption has been based on Rodahl 



and Kasemo's argument 541 ] . which links k to the ratio of force to velocity for low Reynolds' number 
drag on small bodies (e.g. Stokes drag), in which case k is a real constant. This argument ignores 
the fact that important factors in the determination of Stokes drag include the contribution of pressure 
surrounding the particle, the requirement for steady non-oscillatory flow and use of the NSBC. We note 
that friction and drag parameters are usually real because they are not used to describe oscillatory flow. 
In conjunction with this point, we note that the term 'coefficient of friction', which has frequently been 
used to describe k in slip studies, is distinct from the usual parameter of that name, which relates the 
shear force to the normal force in the case of solid-solid sliding friction. Secondly, an argument for phase 



restriction 33J] has drawn on the correlation between bond strength and bond length, which should 
determine magnitude and phase of slip respectively. Such an argument does not explain why the phase 
should take one particular value, or why such a correlation between magnitude and phase should not be 
experimentally investigated. 

Non-zero phase is not significant in previous work relating to lower-frequency oscillations in AFM 



or SFA devices 



231 ] . in which 5 ~100 fxm for 39 Hz oscillations in water, as opposed to ~500 nm for 



551 ] are adequately described by a 



a 1 MHz oscillator. Similarly, periodically patterned microchannels 
real-valued slip length when the ratio of the flow velocity to the pattern wavelength does not exceed 
~1 kHz, corresponding to patterns of extremely short spatial wavelength. 
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B. Slip measurement using the composite torsional ultrasonic oscillator 

1. Theory applied to the torsional oscillator 

The previous discussion of slip is now applied to the composite torsional ultrasonic oscillator (Fig. 
The d amp ing torque when the specimen rod is submerged in a Newtonian fluid was derived by Robinson 
et al. 46J] under the NSBC, using an equivalent circuit analysis. In that derivation, the shear damping 
stress on a flat surface is integrated over the surface of the specimen rod. Appendix A, in this paper, 
gives the same derivation with slip included and Appendix B discusses an assumption that is used. 
Similar derivations have considered dampin g o f other torsional oscillators by viscoelastic fluids under 



the NSBC, yielding similar functional forms 



41 



43, 



561 ] . For the oscillator of specimen rod radius a » 5 



described by Fig. [H the total torque for h(= 1 or 2) immersed flat surfaces is 

T (t) = a 2 ^A + 7^™^ ^ L 'T]a u s,m,ax (t) , (15) 

where the subscript 'max' refers to the value at the antinode of the standing wave, which has wavelength 
A. 

When a specimen rod is immersed, the resistance and inductance added to the equivalent circuit are 
R v and L n respectively, so 

T (t) = (R v + iuL ri ) Us > max {t) . (16) 

a 

Equation \W\ is directly analogous to Eq. [TJ] for this specific configuration if we substitute notations 
R = R v and X = ujL v , consistent with the resistive and inductive descriptions of these parameters. 
Comparing Eqs. [15] and [16] and introducing a constant C, we find 



1 , A V2V , /7T 
l„ — a A H — nira — — a cos — 
2 15 V4 



(7T 
4 ~*~ ^ c 



(17) 



L„ = — a sin I — 

u) V4 
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The equivalent circuit reduces to a series LCR circuit [451 ] . L v is therefore related to the change in 
period upon immersion and the resonant period r, using 



L, 



T '1 

LitrO. • 



(18) 



The inverse quality factor of the oscillation Q 1 is given by the total resistance of the composite oscillator 
Rtot and the total inductance L to f- 



Q 



-i 



R 



tot 



(19) 



tot 



Rri is calculated by assuming that Ltot is dominated by the mechanical inductance (moment of inertia) 



of the rods 



4a], so that (L r) « L 



tot) 



R r] — LdLtotQr) 1 - 



(20) 



Q rj 1 is the change in Q 1 upon immersion. When piezoelectric elements are incorporated, the equivalent 
circuit for the torsional oscillator gives the inverse quality factor, 



(21) 



where Vd and V g are voltages applied to the drive and gauge crystals respectively, and K is a calibrated 
constant for a particular oscillator arrangement [45(. The quality factor is also given by 



Q- 1 



V37 



(22) 



where r is the resonant period and t\ and t 2 are measurements of the period when Vd is double the 
resonance value with V g held constant. By plotting the right hand side of Eq. 1221 against the ratio Vd/V g , 
Eq. [21] can be checked and the value of K determined. K is independent of specimen rod damping, but 
it is dependent on the rigid moment of inertia of the specimen rod, so changes if mass is added to the 
specimen rod. 
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We now consider comparison of a base case (subscript 0) in which the NSBC holds (a = 1, cp a = 0), 
with a case in which slip is present (subscript 1). In practice, it is difficult to determine whether there is 
slip in the base case, so a can be used to measure the relative amplitude and phase of slip between two 
cases investigated. C remains constant between measurements as long as a, A, rj, p and uj are constant, 
so using Eqs. [T71 [18land[20l we define AQ _1 and At as 



Q 
Q 



ILL 

-i 

rfi 



R-ql 
R-qO 



V2\a\ cos + 4>a \ , 



(23) 



and 



At = — = = V2\a\ sin (~ A + 4> a ) . 
T v0 L v0 V4 



Rearranging Eqs. 1231 and 1241 a is calculated directly from measurement: 



w = ,toI + (^!, (25) 



and 



At \ ir . , 

arCtan [aQ^J ~ 4- (26) 



Equations [25] and [26] are plotted in Fig. [5] Figure 5(a) shows that \a\ increases with both AQ~ l and 
At, indicating that if mass is decoupled from the surface, the amplitude of fluid oscillation will decrease 
unless more power is applied. Figure |5(b)| shows that to retain the NSBC value of phase (4> a = 0), 
AQ -1 must decrease if At decreases. Therefore, an oscillation generating relatively low power will 
tend to produce a value of <p a > (fluid oscillation leading) unless there is some decoupling of mass 
from the surface, and therefore reduced inertia. Similarly, if more power is generated, the fluid will lag 
unless balanced by greater inertia. Considered together, Fig. 5(a) and Fig. |5(b)| show that a non-zero 



measurement of either AQ 1 or At necessarily indicates a departure from [a| = 1, or a = or both. 
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2. Viscoelasticity 

The introduction of a slip length is mathematically equivalent to considering a viscoelastic fluid. By 
inspection of Eq. [TTJ and remembering that 5<x^/rj, we find that instead of using a, we can introduce 
a complex viscosity 77* to replace the Newtonian viscosity rj such that a^Jrj = y/rj*. Due to coupling 
to the decaying wave in the fluid, we have aoc^/rj^cc^J where G* = G' + iG" is the complex shear 
modulus. We find that the storage modulus G' and the loss modulus G" are given by combinations of 



the inertial and dissipative terms in the oscillation 



pG' = (R 2 - X 2 ), (27) 

and 

pG" = 2RX. (28) 

It is therefore misleading to use the viscoelastic description of 'storage' and 'loss' moduli for the current 
situation. It is also misleading to describe a Newtonian fluid as viscoelastic until there is experimental 
evidence for viscoelastic behaviour rather than some other cause of slip. Some previous QCM work based 



36. 



57|, 



581 ] has considered slip in conjunction with 



on fully phase-dependent equivalent circuit analyses 
a viscoelastic fluid and layering, thereby confusing possible Newtonian slip with other effects; for further 
effects of layering see the Discussion (below). Such an approach may be motivated by layering of a QCM 



34] 



electrode substrate, that is not present when using a torsional oscillator. One previous investigation 
has noted the equivalence of viscoelasticity and slip in interpretation of slip measurements. 

3. EXPERIMENT 

A. Experimental details 

To demonstrate investigation of slip using a torsional oscillator, preliminary experiments have been 
carried out using a 40 kHz oscillator, as described by Fig. SJ The fused quartz specimen rod was 
uniformly ground using AI2O3 grits to remove surface features of magnitude >5 pm, and polished using 
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cerium oxide powder. Polishing debris was removed by sonication in acetone, isopropanol, 10% nitric acid 
and ethanol. Any remaining organic material was removed using 'piranha' solution (70 vol% sulphuric 
acid, 30 vol% hydrogen peroxide). The rod was then thoroughly washed in deionised water (0.2 fim 
filtered, 18.2 MO cm) and lightly flame-polished to minimize surface roughness. The 'slip' case was 
configured by coating the specimen rod with a hydrophobic fluoropolymer surfactant (RS Components), 
evenly dispensed from an aerosol can. The base case was retrieved by dissolving the surfactant in 
isopropanol prior to re-cleaning the rod using 'piranha' solution and washing with deionised water and 
ethanol. 

A pre-optically polished, flame-polished piece of fused quartz has RMS roughness < 1 nm and peak 
to trough roughness < 3 nm, as measured over several widely-distributed areas of size 5 ^m 2 using 
AFM. Previously 



471 ]. Scanning Electron Microscopy (SEM) revealed peak to trough roughness of 2— 
4 fj,m on rods left unpolished following grinding using AI2O3, with no significant alteration of the texture 
following the addition of the surfactant. However, a flame-polished surface coated with the surfactant 
has additional texture characterised by features of spacing < 5 fim with RMS roughness 24 nm and 
peak to trough roughness < 300 nm, observed using optical microscopy and AFM. The advanced static 
contact angle of deionised water on fused quartz is 44 + 2°. The advanced contact angle on the same 
surface coated with the surfactant is 118 + 1°. The texture of the coating causes some hysteresis, and the 
receded contact angle is 113 + 1°. Contact angles were measured using the sessile drop method (CAM 
200, KSV Instruments). 

Prior to each immersion run, the specimen rod was washed in deionised water and dried in air. For 
the base case, the uncoated rod was also rinsed with ethanol prior to drying. Deionised water was used 
as the test fluid in all experiments. The oscillator was contained in an insulated bell jar standing on 
a plate maintained at 303 + 1 K and situated on a large mechanically damped pad usually used for a 
Transmission Electron Microscope. 
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B. Measurement and calibration 

Measurements from recent work are presented alongside previously reported experiments [^tJ in Table [J 
and explained in this section. The data points for AQ _1 and At are each the result of several immersion 
events in each of the base and slip cases. The specimen rod was immersed by slowly raising a filled vessel 
to the rod. A value of Q -1 or for any particular immersion run is the difference between the dry value 
and the value when the water was in contact with the node of the standing wave, adjusted to account for 
any background drift. The fluid meniscus was aligned with the vertical mid-point of a rubber o-ring that 
was fitted half way along the specimen rod. The o-ring, which fixes the position of the node, was located 
at a point where surface motion is zero or minimised, and is considered rigid during calibration. Each 
immersion was carried out after the oscillator had been left for at least 90 minutes, so that instrumental 
drift per second was < 5e -5 for Vd/V g and < le -13 s for r, as measured over 10 minute intervals. 

The standard error of the mean measurement was combined with error ascribed to drift to give the 
total experimental error in Q~ l and in each case. Uncertainty is dominated by variation between 
different immersion events, probably due to the extreme thermal and mechanical sensitivity of the 
instrument. The lesser contribution due to instrumental drift is accounted for by assuming that the 
mean drift correction is a random error. Earlier experiments have established that uncertainty due 
to surface tension effects and location of the node is small relative to the standard error in the mean 
value, because very little damping occurs at or near the node. 

In order to calculate Q^ 1 from measured voltages, the half width of the oscillation was measured 
directly following each immersion. Equations [2H and [221 were then used to calculate a precise value of 
K. It has been shown \^\ that a plot having gradient K is typically linear (R 2 > 0.9999 for all fits) 
and independent of temperature, the immersion level and the presence of surfactant. The two former 
variables are held constant during immersions, but the slip case experiments were calibrated separately 
from the base case in the most recent experiments. The value of K is slightly lower for the slip case, 
as expected from the first principles analysis 



461 ] in which K is inversely proportional to the oscillator's 
moment of inertia. The maximum error in K, conservatively estimated as the difference between the 
slip case and the base case, is not significant in comparison with the experimental error of a. K changes 



10 



when the torsional oscillator is mechanically altered, but was constant for each individual shear rate 
presented in FigO 



For the NSBC, it is expected that Q" 1 = 2t v /t 46|. This relation was followed to within 5% for 
the majority of immersions. Even accounting for the possibility of slip, immersion events with >10% 
deviation clearly indicate a significant instability in the oscillation and such data were discarded. 

Experiments have been performed at shear rate amplitudes \i r $\ between 230 and 6840 s _1 , corre- 
sponding to linear displacement amplitudes of 3.2 and 96 nm at the antinode. The linear amplitude 
of the oscillation is equal to_V25\i r g\/io. The shear rate at the surface of the rod is calculated as de- 



scribed by Robinson et al. 45J]. For this calculation, rod diameters and masses were measured using a 
micrometer screw gauge and a Mettler AE 200 mass balance, respectively. For each strain rate, intrinsic 
gauge capacitance C m was measured to better than +10% by finding the gradient of straight line plots 



45(]. 



(R 2 > 0.999) of the fractional change in resonant period as a function of total gauge capacitance 
The transformer ratio linking the electronic and mechanical parts of the equivalent circuit was calcu- 
lated using K. The alternative calculation (from first principles) incorporates large error relating to the 



electrode separation angle 2$ (Fig. 01 46J]). The final uncertainty in shear rate amplitude is ~ 15%. 
However, the uncertainty in each value relative to the others depends only on V g , which is precise to 
better than 1%. 

C. Results 



The results for the five experiments in Table Hare plotted as slip parameters in Fig. [6j Figures 6(a) and 
|6(b)| appear similar because, for all five data points, mean values of (AQ _1 — l) and (At — 1) have the 
same sig n. All data lie within experimental uncertainty of the NSBC. For the three data points from 
Ref. [47], the consistency of mean values (|a| < 1, <p a < 0) previously led to the conclusion that finite 
slip had likely been observed. Although the newer results are not directly comparable due to differences 
in surface preparation, the mean values are now scattered around |a| = 1, <j) a = 0, suggesting that the 
NSBC is observed with a high degree of accuracy. The most recent experiments have covered a greater 
range of strain rates and afforded greater precision due to improved stability of the oscillation and better 



17 



thermal and mechanical isolation. The corresponding mean values of |6| are now less than 30 nm, and 
the upper limits bounded by experimental uncertainty are less than 60 nm. For each data point, the 
lower error bar of |6| is at zero because the NSBC falls within experimental uncertainty. Because of the 
phase discontinuity at the NSBC (Fig. [3|), any value of 4>b is possible for near-NSBC combinations of 
the measured variables. 

Previously, QCM measurements 



33f | have obtained values of b of similar precision, although in that 
case finite slip lengths of magnitudes between 10 and 80 nm were observed on bare as well as a range 
of hydrophobic surfaces. As discussed previously, the analysis used to derive these values contains 
inconsistencies. When considering other results in the literature, we emphasize that comparisons are 
drawn from data sets taken from measurements of very different flows (see Discussion) and that the 
current data are observed over a well-defined, wide range of shear rates. It is therefore remarkable 
that the data is consistent with emerging trends from high-precision AFM studies. Namely, slip at a 
hydrophobic smooth surface occurs on scales < b = 50 nm, if at all; and surface roughness such as that 
observed for the slip case in the current work may play a part in suppressing any slip at a hydrophobic 
surface. 



4. DISCUSSION 

A. Interpretation of slip measurements 

In order to discuss the interpretation of experimental slip measurements, we first clarify the terminology 
used in the discussion. Lauga et al.'s definitions Q| of 'indirect' measurement, 'apparent' slip and 
'effective' slip are used, although we extend the definition of effective slip to include temporal as well as 



.Q. 



We define it here as 



spatial averaging. 'True' slip is an oft-used term not defined by Lauga et al. 
the case when the mean tangential velocity component of fluid is different from that of a solid surface 
immediately in contact with it. 

Experimental methods are predominantly indirect at present. Even the most 'direct' optical methods 



of slip measurement have been labelled indirect 



23J , in the sense that the velocity of tracers rather than 
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the fluid itself is measured. In the case of the oscillator, the measured damping force is not necessarily 
consistent with the idealised analysis presented in the Theory section. In any surface force measurement, 
the flow profile is inferred rather than explicitly determined. Under these conditions, slip parameters b 
and a only describe the slip dynamics in full when rate- independent, true slip occurs and the fluid is 
entirely Newtonian. 

Caution is required when analytically linking postulated mechanisms to apparent or effective slip 
measurements. For example, it often appears convenient to assume the presence of layering at or near 



the slip interface 
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3a 



36 



38 



39. 
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50 



58|| . Typically, layers of unknown thickness and interfacial 
properties only introduce more modelling parameters. As discussed previously, introduced viscoelas- 
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351 ]. multiple viscoelastic layers 



ticity is indistinguishable from slip. In two studies employing a 
with multiple boundary conditions were modelled and several parameters were fitted. The resulting 
magnitudes of a (ranging from ~ 1.8 to ~ 5) were larger than reasonable for a Newtonian fluid. At 
present, mechanisms of apparent slip are poorly characterized, so empirical measurement of slip param- 
eters remains the most important step. Newtonian fluids should be considered Newtonian right up to 
the interface to avoid introducing error. 



B. Comparison with other experimental techniques 

SFA and AFM techniques are currently the most well developed and widely-employed methods of mea- 



suring surface slip, and can reach resolution of better than +5 nm for b 
as the Fluorescent Recovery After Photobleaching (FRAP) 
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24|. Optical methods such 



251 ] and Particle Image Velocimetry (PIV) 



591 ] techniques allow resolution of +100 nm, limited by diffraction and tracer particle size. Recently, 



a nano-PIV system has reportedly resolved true slip lengths of less than 100 nm with uncertainty less 
than +20 nm Several groups have measured slip using capillary techniques [l|, Q], but accuracy 

has been limited by capillary fabrication issues, the difficulty of controlling surface chemistry within a 
capillary and (in some cases) the effects of electro-osmotic flow. A recently-developed technique derives 



slip from measurement of thermal motion of confined colloidal particles 



. A slip length of 18 + 5 nm 



has been reported at near-zero shear rate. Torsional oscillator experiments presented in the current 
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work confirmed the NSBC to within +0.03 for \a\, +1.5° for 4> a and to within 60 nm of |6| = 0. We 
stress that there are no significant limitations on the precision of the oscillator other than control of 
atmospheric conditions, vibration and electrical noise. Further, the technique offers advantages and 
points of difference when compared with other techniques. 

A particular advantage of the torsional oscillator is characterization of experimental shear rate. The 
shear rate amplitude is easily varied by adjusting the applied gauge voltage and can at least cover the 
range 10 2 to 10 4 s- 1 . A major reason for use of torsional oscillators in rheological experiments is the 



large and difficult-to-reproduce range of accessible shear rates 



4& 



441 ]. In contrast, the variable vel ocity 



34J. In the AFM and SFA techniques 



201, 



distribution across a QCM causes analytic complications 
the shear rate varies with respect to surface separation and distance from an axis passing through the 
centres of the approaching curved surfaces. The shear rate is maximized away from this axis and also 



3C]. 



depends on the approach or withdrawal rate of the two spheres 

The second major advantage of the torsional oscillator is versatility. Surface treatments can easily 
be applied to the specimen rod, either while remote from the oscillator or attached to it. If the specimen 
rod is detached, the response of the oscillator is easily recalibrated. Electronic voltage and period mea- 
surements are captured near-instantaneously, so time-dependent trends can be measured. The oscillator 
is perfectly suited to measurement of in-situ changes of surface chemistry or other parameters, and 
experimental precision would be significantly enhanced for such experiments. The oscillator can readily 
be used at variable temperatures and even variable pressures. Although it is best suited to differential 
measurements (similarly to optical techniques), the zero-slip response can be determined by careful cal- 



ibration of all the parameters described in the original equivalent circuit analysis 
use of QCMs and other torsional oscillator designs 
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43 



45 



461 ] . In contrast, 



561 ] for slip measurement presents difficulties 



with regard to surface modification and incorporating piezoelectric equations into an equivalent circuit 
analysis. 

An important distinction between damping-oscillation methods and other techniques is that the 
shear damping force is measured. SFA and AFM methods measure drainage forces caused by pressure 
in the flow around curved surfaces; viscous forces are not significant. Another difference is that, when 
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using the torsional oscillator, slip occurs over a relatively large fluid-solid interface, so measurement of 
molecular-scale effects is limited by surface homogeneity. However, measurement of effective slip will give 
more accuracy and precision than a value derived from a small system. It is not true that measurements 
of effective slip must take place on a scale L ~ b in order to have a measurable impact Q], because any 
fundamental slip mechanism should be consistent over the whole interface. Regardless, we note that the 
linear displacement amplitude for the torsional oscillator is of comparable magnitude to the slip length, 
so molecular-scale interactions are certainly probed using this device. 

5. CONCLUSION 

In this paper, we have drawn particular attention to the occurrence of slip at an oscillating interface. The 
ultrasonic composite torsional oscillator is a versatile instrument for measurement of viscous damping 
forces, and therefore slip. In order to fully characterize this device for slip measurements, we have 
revisited fundamental aspects of the slip analysis. The analysis has then been applied at a generalised 
oscillating surface, followed by the surface of the oscillator's specimen rod in particular. Preliminary 
experiments at shear rates between 230 and 6800 s _1 have probed slip at a hydrophobic surface of RMS 
roughness 24 nm, obtaining resolution better than +0.03 for \a\, +1.5° for <p a and within 60 nm of 
\b\ = 0. 

Navier's slip parameter, while retaining the units of length, is best suited to interpretation as a 
surface force effect. Measurement of surface forces has played an important role in developing the 
understanding of slip to date, while applications are likely to utilise the lubricating effect at a slipping 
surface. The actual velocity discontinuity at an oscillating interface, as encapsulated by the parameter a 
for oscillating surfaces, is being directly observed with increasing precision using optical methods. The 
link between a velocity discontinuity and a surface force necessarily utilizes the spatial gradient of an 
analytical flow profile, unless both types of experiment can be performed on the same system. 

Ongoing experiments will aim to significantly develop the stock of published slip measurements at 
oscillating surfaces. Particular experimental topics will be the influence of shear rate, hydrophobicity and 
surface texture, while the development and in-situ characterisation of surfaces or fluids with switchable 
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properties will also be a priority. 
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APPENDIX A: DERIVATION OF SLIP-ADJUSTED TORQUE 
APPLIED TO A SPECIMEN ROD 

In this section, the damping torque on a composite torsional ultrasonic oscillator rod, as illustrated in 
Fig. 21 is calculated. We follow Robinson and Smedley [46f] and incorporate the case in which there is 



finite slip, as introduced in the Theory section. Using the cylindrical polar co-ordinate system defined 
in Fig. U the velocity in the azimuthal (6) direction is given by 

u e (r, z, t) = u s , max (t) cos f-y- J , (Al) 

where the subscript 'max' refers to the value at r = a and at the antinode of the standing wave, which 
has wavelength A. 

There is a velocity gradient in the fluid in both the r- and z-directions. Viscous forces in the fluid 
exert azimuthal shear stresses on the rod, which dampen the torsional oscillation. Forces due to the z- 
gradient are neglected (Appendix B). The total damping torque is derived by integrating the tangential 
stress (Eq. [TT]) along the immersed, curved length of the rod I. It is assumed that the fluid surface 
coincides with a node and the free end is an antinode, so that sin (^) = 1. The torque is 

Trod (t) = <?er (t) 2-Ka 2 dz 

2 ^-7ra 2 r]au s>max (t) J cos^^^^cfe ( A2 ) 



6 

a 2 \r]au S) max (*) • 



f2% 2 
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The torque on the oscillating flat end of a rod is derived similarly. The surface velocity at radius r 
from the cylindrical axis is given by 



u g (r, z = 0, t) = -u s max (t) , 
a 



(A3) 



and the torque is then integrated over the surface of the free end, so that 



Tend 



(*) = f ae» (t) 
Jo 



27rr 2 dr 



2V2i 

5a 
\/2i7ra 3 



(t) r 3 dr 
Jo 



(A4) 



6 2 

Neglecting edge effects between the flat end and the curved surfaces, the total torque for h immersed 
flat surfaces is 



T (t) = a 2 (A + ^hira)—z-r]au Stmax (t) 



(A5) 



APPENDIX B: VISCOUS GRADIENT IN THE z-DIRECTION 

The velocity flow field around the torsional rod has gradients in both the r- (radial) and z- (axial) 
directions. The azimuthal stress due to these gradients are given by 



a B r = rj- 



du e (r, z, t) 



(T 



(Bl) 



-u e (a,z,t) 



and 



0"0 2 (t) = V 



du e (r, z, t) 



2tt 

"a 



dz 



ug (a,z,t) . 



(B2) 



For a 40 kHz torsional oscillator in water, the ratio of radial to axial stresses is ~ W 4 
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FIG. 1. The slip lengths b and b\ defined in the text, with the shaded area representing solid material. 
Navier's slip length b (Eq. [3]) is only equal to the variable b\ (Eq. [7]) when the velocity profile is linear. 
The labels x and z denote Cartesian coordinate axes. 
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FIG. 2. Magnitude (Fig. 2(a)) and phase (Fig. 2(b)) of u»/u s as a function of the distance from an 



oscillating surface for a Newtonian fluid, calculated using Eq. [9] and plotted for the NSBC as well as 
various slip boundary conditions. The NSBC phase is discontinuous at the interface, taking a value of 
zero at z = even though the limit as z/5 approaches zero is 45°. 
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FIG. 3. Magnitude (Fig. 3(a) ) and phase (Fig. |3(b)[ ) of a relative to b for a Newtonian fluid, calculated 
using Eq. fTUl Equations [3] and [9] can be used to show that <f>b takes a value of —45° rather than 135° 
when (p a = 0. A value of (f>f, = 135° is consistent with a negative value of b and NSBC phase dependence. 
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45|, 



461 ] . Figure 4(a) shows a piezoelectric x-cut 



FIG. 4. The composite torsional ultrasonic oscillator 
a-quartz rod, which takes dimensions of height H x 48mm and of radius a ss 5 mm for the rods of 
resonant frequency 40 kHz used in experiments presented here. Four conducting electrodes, separated 
from each other by an angle 2$, are coated on to the rod surface and a voltage V is applied in order 
to excite torsional oscillations. In the composite oscillator (Fig. |4(b)[ ) two identical piezoelectric rods 
are joined together with a co-resonant specimen rod. A standing torsional wave of wavelength A is set 
up with antinodes located at free ends and joints. Typically the length of each rod is A/2, as in the 
figure. The apparatus is held at the stationary nodes, where the electrode connections are also located. 
Viscous fluid damping is measured by raising the test fluid on an adjustable stage so that the specimen 
rod is immersed from the free end along a length of the rod I. The z- and r-axes of the cylindrical polar 
coordinate system used in the text are as shown. The voltages across the drive (V^) and gauge (V g ) 
crystals are measured along with the resonant period r. 
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FIG. 5. Magnitude (Fig. 5(a) ) and phase (Fig. |5(b)[ ) of a for a Newtonian fluid, calculated using Eqs.l25l 
and[26l Each line is plotted at a constant value of either AQ -1 or At, with the other measured quantity 



taking the value determined by the x-axis. In Eq. [25] and Fig. 5(a) the two measured variables are 
equivalent. 
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FIG. 6. Magnitude (Fig. 6(a) ) and phase (Fig. |6(b)[ ) of a and magnitude of b (Fig. 6(c) ) for the 



experiments in Table [H Va 
those presented in Ref. 



ues of |6| have been updated slightly (< 20 nm difference in all values) from 



471 ] so that drift is consistently handled. Error bars for the calculated shear 



rates are omitted for clarity. 
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TABLE I. Details of experiments previously presented ( 43], labelled 'A') and new experiments (labelled 
'B'). 



Experiment 


£r6> 


Linear Amplitude 




K 


AQ 


-i 


At 




(10 3 S" 1 ) 


(nm) 




(io- 5 ) 








A-l 


0.43 ± 0.06 


6.1 ±0.9 




0.93 


0.99 ± 


0.01 


0.95 ±0.05 


A-2 


1.4 + 0.2 


19 ±3 




0.87 


0.99 ± 


0.01 


0.94 ±0.07 


A-3 


3.6 ±0.5 


51 ±7 




0.36 


0.99 ± 


0.02 


0.96 ±0.10 


B-l 


0.23 ±0.03 


3.2 ±0.5 


(base) 


3.60 


1.00 ± 


0.01 


1.02 ±0.03 








(slip) 


3.59 








B-2 


6.8 ± 1.0 


96 ± 14 


(base) 


3.57 


1.00 ± 


0.04 


1.00 ±0.01 








(slip) 


3.57 
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